Efficient Collection of Single Photons Emitted from a Trapped Ion into a Single Mode 
Fiber for Scalable Quantum Information Processing 
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Interference and coincidence detection of two pliotons emitted by two remote ions can lead to an 
entangled state which is a critical resource for scalable quantum information processing. Currently, 
the success probabilities of experimental realizations of this protocol are mainly limited by low 
coupling efficiency of a photon emitted by an ion into a single mode fiber. Here, we consider two 
strategies to enhance the collection probability of a photon emitted from a trapped Yb+ ion, using 
analytic methods that can be easily applied to other types of ion or neutral atoms. Our analysis 
shows that we can achieve fiber coupling efficiency of over 30 % with an optical cavity made of a flat 
fiber tip and a spherical mirror. We also investigate ways to increase the fiber coupling efficiency 
using high numerical aperture optics, and show that collection probability of over 15 % is possible 
with proper control of aberration. 

PACS numbers: 32.80.Qk,42.50.Pq,42.15.Eq 
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INTRODUCTION 



Entanglement is a unique feature of quantum physics 
and it has found various applications ranging from the 
experimental verification of the foundation of quantum 
physics [l| to many quantum information protocols suj;h 
as quantum teleportation J3| , quantum cryptography '3] , 
quantum computation Q , and long distance quantum 
communication [6]. Most of these applications require 
a large number of entangled states distributed over long 
distances, which can be generated using one of two ap- 
proaches. In the first approach, one can locally gener- 
ate a pair of entangled photons 7], distribute them to 
remote sites, and map these flying qubits into station- 
ary (memory) qubits Q. To efficiently transfer a flying 
qubit into a stationary qubit, the spectral and spatial 
modes of the flying qubit have to be matched to the time- 
inverted emission properties of the memory qubit, which 
is the most significant challenge in this approach. In the 
second approach, one first generates flying qubits that 
are entangled with a local stationary qubit at each site, 
sends these photons to a common location, and then in- 
terferes these two photons using a 50:50 beam splitter 
d, With 50 % chance, these photon pairs are pro- 
jected into one of two distinguishable Bell states, and 
once one of these two Bell states is detected, the two sta- 
tionary qubits at the remote sites are also projected to 
a maximally entangled Bell state The advantage of 
this approach is that the spectral properties of the two 
photons are naturally matched if one uses the same type 
of memory qubits. In order to achieve interference of the 
two photons with high visibility, the photons have to be 
indistinguishable. Therefore, we need to collect emitted 
photons into a single spatial mode, such as a single mode 
fiber. Several experiments demonstrated establishment 
of entanglement between remote trapped ions using this 
approach [T2| - |l^ |. While these experiments realized high 
fidelity entangled states, the success probability of the 



probabilistic entanglement process is very small. The 
success probability is mainly limited by the small photon 
collection probability of about 0.4 % given by the small 
numerical aperture (NA) of the collection optics and an 
insufficient overlap between the spatial mode of the single 
mode fiber and the emission pattern of an ion. Because 
this protocol requires one to collect photons from both 
sites and observe a coincidence, the overall success prob- 
ability scales quadratically with the single fiber coupling 
efficiency, and was limited to only about 2 • 10^^ (see Ta- 
blc[II]in Sec. IIVI) (l^ . To apply this system for large scale 
quantum information processing or long distance quan- 
tum communication, it is critical to dramatically enhance 
the collection probability of photons into a single mode 
fiber. 

In this paper, we show systematic approaches to in- 
crease the fiber coupling efficiency of a single photon 
while preserving ion-photon entanglement. We will an- 
alyze the design of an optimized optical cavity that can 
be integrated with a surface ion trap, and evaluate vari- 
ous ways to achieve high NA collection optics compatible 
with surface traps to enhance light collection into a single 
mode fiber. 

Our analysis mainly focuses on ^^^Yb''' whose sim- 
plified level scheme is shown in Figure [I] (a) pj, but 
the conclusions can be applied to other atomic systems. 
There are several ways to entangle a photon's degree 
of freedom with internal electronic states of a trapped 
171YI3+ [qji jf Zeeman levels are denoted by \F,mp), 
where F is the total angular momentum and ttif is its 
projection onto the quantization axis, one can excite the 
Yb ion to '^Pi/2\^,0) and collect the photon from the 
(T* transition as shown in Figure [1] (b). The polar- 
ization of the photon will then be entangled with the 
two ground state Zeeman levels ^5*1/211, ±1) [l^. Al- 
ternatively, one can collect a photon generated by a- 
transitions from ^Pi/2|1;~1) as shown in Figured] (c) 
or by TT-transitions from a superposition of ^Pi/2|1,0) 
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FIG. 1. Ion-photon entanglement schemes in ^^^Yb"''. (a) 
Simplified level scheme of ^^^Yb^ with dipole-allowed tran- 
sitions. Zeeman levels are denoted by \F,mF), where F is 
the total angular momentum including nuclear spin and mp 
is its projection onto the quantization axis. ^S'i/2|1,0) and 
■^5i/2|0,0) are first-order magnetic field-insensitive hyperfine 
levels (clock states), (b) When the ion decays fi'om ^Pi/2|1, 0) 
to one of the Zeeman levels in ^Si/2|1, ±1), the polarization 
of the emitted photon gets entangled with two Zeeman levels 
. (c) The frequency of the emitted photon and the atomic 
clock states are entangled after a decay from ^Pi/2|1, —1) 
(d) Due to the selection rules, excitation of ^^^Yb"*" prepared 
in a superposition of the clock states in ^S'1/2 generates a co- 
herent superposition of the excited states (two states 
with TTiF = 0). Once excited state decays back to the original 
state, the frequency of the emitted photon is entangled with 
clock states In ^'^^Yb+ system, ^51/210,0) is generally 

used as an initial state due to its high preparation fidelity. 
Therefore (c) and (d) scheme require additional microwave or 
optical Raman transitions denoted by thick arrows. 



and ^^"1/210,0) to the atomic clock states composed of 
25i/2|l,0) and 25i/2|0,0) as shown in Figure[T] (d). The 
frequency of the photon is then entangled with the clock 
state qubits. The spatial emission pattern of photons 
from and 7r-transitions and their spectral proper- 
ties together with the ability to filter out photons from 
unwanted transitions will dictate strategies that can be 
adopted to enhance the photon collection efficiency. For 
remote entanglement generation, we will analyze the ap- 
proach using an optical cavity in Section and an 
approach using high NA optics (either a refiector or a 
lens) in Section IIIII Section IIVI estimates the entangle- 
ment generation rate we can achieve based on these ap- 
proaches, and conclusions are presented in Section Ivl 



II. CAVITY ENHANCED PHOTON 
COLLECTION 

A small high finesse cavity can be used to collect light 
from a single atom by induced emission into the cavity 
mode. Although efficient light collection using an optical 
cavity has been demonstrated for neutral atoms [3) 113 > 
additional challenges have to be met for trapped ions: 
readout transitions of trapped ions are typically in the 
ultraviolet (UV) and the ion trap has to fit in a small 
cavity with UV-dielectric coatings. In addition, trapped 



FIG. 2. Geometry of an optical cavity integrated with surface 
trap. Optical cavity is composed of a small spherical mirror 
and a fiber in a ferrule whose tip is coated with high refiecting 
dielectric layers. I is cavity length, wq is beam waist at the 
fiber tip, and rton is beam radius at the trapping position 
of the ion. Surface trap can be patterned with Au on the 
fiber/ferrule surface whose diameter is >2 mm. 



ions have to be shielded from the charges which can ac- 
cumulate on dielectric surfaces p^ . 

The collection probability of a single photon from an 
excited atom into a single mode fiber is dictated by sev- 
eral factors. First, the excitation has to leave the coupled 
atom-cavity system via the cavity decay channel rather 
than spontaneous emission, characterized by the photon 
extraction probability PcavUy ■ While some photons leave 
the cavity through the desired output channel, other pho- 
tons are either lost due to scattering and absorption in 
the mirrors or transmitted through the undesired mir- 
ror (s). We define rt as the fraction of the desired trans- 
mission over the total cavity decay. Finally, we consider 
overlap between the mode of the cavity and that of the 
single mode fiber, which will be called mode matching 
factor e. The overall probability P fiber of extracting the 
photon through the single mode fiber is then given by 

P fiber — ^ ' '^t ' Pcavity • (1) 

In this section, we will investigate how to maximize each 
factor within the physical constraints. We will mainly 
consider a cavity geometry shown in Figure [2] [l9| and 
optimize system parameters to maximize the collection 
efficiency under the given constraints. This discussion 
will also identify components that deserve further inves- 
tigation to improve the performance of the system in the 
future, and provide a metric to compare different cavity 
systems. 

When an atomic system is coupled to a cavity mode, a 
single photon can be generated either by fast excitation 
with resonant 7r-pulse followed by spontaneous emission 
p^ - ffl [20| , or by stimulated Raman adiabatic passage 
(STIRAP) process [IB, In this section, we will also 
compare these two methods in terms of collection prob- 
ability Pfiber and repetition rate. 
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A. Fast excitation with resonant vr-pulse: Design 
principles 

In this subsection, we consider the case where the ion 
is excited by an ultrafast pulse whose duration is short 
compared to any time constants of the system, such that 
the initial state is an excited state of the ion. When an 
excited ion is placed inside an optical cavity, the sponta- 
neous emission rate into the cavity mode can be enhanced 
by the Purcell factor, 2C + 1 [IH, where C = g'^/2K-f 
is the cooperativity, g is the coherent ion-field coupling 
constant, k, is the cavity-field decay rate, and 27 is the 
spontaneous emission rate. 

The parameter is inversely proportional to the cav- 
ity mode volume and k = tt • FSR/T, where FSR and 
T are the free spectral range and the finesse of the cav- 
ity, respectively. Therefore, we need a high finesse cav- 
ity with small mode volume to achieve large coopera- 
tivity. Small mode volume cavities can be realized by 
many different approaches such as microdisk , micro- 
sphere [2^, microtoroid [13], and photonic crystals [25j . 
Integration with ion traps requires a cavity with empty 
mode volume, so we will consider only traditional open 
cavities formed by two (curved) mirrors. Cavity mode 
shape between two spherical mirrors is Gaussian and the 
mode volume is 7rwg//4, where wq is the beam waist of 
the Gaussian mode function and / is length of the cavity. 
The cooperativity is reduced to C = 3J"-Br- A^/(7r^rfg„), 
where Br is a branching ratio, A is wavelength of transi- 
tion, and rion is the radius of the mode function at ion 
location. In this system, we need high T and small rion 
to maximize C. 

Using these cavity parameters, one can express the 
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where Ic is a characteristic cavity length defined as Ic — 
TTc/(2jJ-) and Vc is a characteristic beam radius defined 
as Tc — ^y6X^J^ ■ Bt/tt^. Eq. ([2]) implies that if either 
I or Tion is much larger than the corresponding charac- 
teristic lengths, the overall Pcavity will decrease rapidly. 
The characteristic cavity length (beam radius) decreases 
(increases) as the finesse of the cavity increases. This 
indicates that for high finesse cavities one should use a 
short cavity but the beam radius can be larger, and the 
design shifts to longer cavities with smaller beam radius 
as the finesse of the cavity decreases. 

Cavity finesse J" is defined by 27rVf — C/C, where 
C = Tf + Te + Lf + Le is the total loss of the cavity, 
Tf and Lf (Tg and Lg) are transmission and passive loss 
through dielectric coating on the fiber tip (the other end 
mirror) in Figure [51 respectively. Passive loss is due to 
scattering loss mainly determined by surface roughness 



[27j and absorption in the dielectric coating. We can 
define Jq = — {Lf + Le)/{Lf + Le), to denote the 

maximum achievable finesse limited by the uncontrollable 
passive loss. Then the actual finesse can be approximated 
to ~ -^0(1 ~ *"«) where rt = Tf /C. Assuming the mode 
matching factor e is constant, one can see that P fiber 
defined in Eq. ([1} is maximized when 



n = no 



l + I/v/(l + 2Co)(l + /Aco) 



(3) 



to a value 



P 



Co 



fiber.max 



I + C0 + l/2l,„ + V(f + 2Co)(f + /Aco)' 

(4) 

where Co = STq ■ Br • X'^ /(TT^rfon) and Ico = t^c/{2jTo) 
are the cooperativity and characteristic length with max- 
imum possible finesse J^o, respectively. From Eq. (|4]), one 
can see that Co should be maximized while I should be 
minimized in order to maximize the fiber coupling effi- 
ciency. 

Finally, fiber collection efficiency is further reduced 
by mismatch between cavity mode {E) and fiber mode 
{G fiber) at the fiber tip characterized by the mode match- 
ing factor 



\J^daE*-Gf, 



her 



^^da\E\^^^da\Gf, 



(5) 



where G fiber and E are the transverse Gaussian mode 
shapes of fiber and cavity, respectively. The optimal wq 
of the fiber cavity should be determined to maximize the 
coupling P fiber to an ion at height /liom which may not 
necessarily match with the fiber mode. If one uses addi- 
tional optics to couple the cavity mode into a fiber mode 
rather than integrate the fiber mode directly underneath 
the flat cavity mirror, the mode matching factor given 
by Eq. ([S]) can be improved substantially, approaching a 
value of 1. 



B. Cavity design example 

The optimum design of the cavity-ion trap system is 
often constrained by the ability to achieve high perfor- 
mance components. In our system, the dominating con- 
straint is the finesse J-q and the radius of curvature (RoC) 
for the curved cavity mirror one can achieve, limited by 
the coating quality and polishing capability. Compared 
to near-infrared wavelengths used in neutral atom exper- 
iments [17} , the passive loss of cavity mirror coating 
is significantly higher in the ultraviolet range relevant for 
most trapped ions. Scattering loss from a rough surface 
IS given by 1 -exp[-(47rcr/A)2] ^ where (7 is RMS sur- 
face roughness. For example, we observed ~ 150 parts 
per million (ppm) loss at 369.5 nm on a super-polished 
substrate whose RMS surface roughness is measured to 
be <0.1 nm, corresponding to estimated scattering loss 
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of <12 ppm. We can attribute most of the loss in this 
case to absorption loss in the coating material in the UV 
wavelength range. Since the RMS surface roughness of 
the fiber tip and the mirror substrate with small RoC is 
generally larger than that of a super-polished mirror sub- 
strate with large RoC, it is reasonable to assume that the 
uncontrolled loss in our mirrors is Lf + L^ wl500 ppm, 
corresponding to ~ 4200. 

Next, we have to pick the radius of the cavity mode 
Tion at the ion location and the cavity length /. In order 
to maximize Cq, one should find a minimum value of 
rion- Since the ion cannot be trapped at the beam waist 
in a fiber tip cavity shown in Fig. ^ , the optimal waist 
Wq is determined by the distance between the ion and 
the waist which is located on the flat fiber tip. When 
this distance is given by /lion, Tion can be minimized if 
uiQ = \Jhion • A/tt. In a surface trap, the minimum ion 
height from the surface is constrained by considerations 
such as the exposed dielectric material |18|, anomalous 
heating [2^ and the scattering of laser light from the 
surface, and /i^on = 50 /xm is a good compromise. In this 
case, the optimal waist is vjq — 2.4 /im, corresponding to 

Tion = 3.4 [IVa. 

For a cavity created between a flat mirror and a con- 
cave mirror, the beam waist of the Gaussian mode is 
located at the flat mirror surface. The length of the cav- 
ity I is related to the Rayleigh range zr = ttwq / A of the 
Gaussian mode as = (RoC — To trap the ion inside 
the cavity volume, we require I > hion = zn and there- 
fore I should be chosen to be very close to the RoC. A 
practical mirror one can polish using conventional super- 
polishing technique has a RoC > 5 mm, so we choose 
I Ri RoC — z^/RoC « 5 mm. The cavity length must 
be aligned to be about 0.5 fim shorter than the RoC of 
the mirror. Laser ablation can create a smooth curved 
surface with RoC between 40 /im and 2 mm with 0.2 nm 
RMS roughness [2^ , but the aperture diameter of a mir- 
ror fabricated in this way is limited to < 100 /um which 
is too small for near-concentric cavity geometry. 

Using these parameters, we obtain Co ~ 4.8 and ^co ~ 
1.8 mm. The transmission of the out coupling mirror is 
determined by setting « 86% according to Eq. ([3]), 
corresponding to Tf « 9,500 ppm. This choice deter- 
mines critical cavity parameters as = 570, C = 0.65, 
Zc=13 mm and = 3.9 /im. Measured beam waist of a 
typical UV flber at 369.5 nm is about 1.5 iim and in Fig- 
ure [3] the change of the Pf function of the cavity 
length I for an ion height hion — 50 /im is shown. The 
mode matching efficiency e is about 82 % if wq is opti- 
mized for maximum Pcavity {wo = 2.4 /J.m), and Pcavity is 
reduced when e is 100 %. Therefore, there is an optimal 
length that maximizes P fiber , near which one can achieve 
more than 30 % of coupling efficiency compared to 0.4 % 

If additional optics is used to couple the cavity mode 
into a single mode fiber (instead of using fiber tip as 
the fiat mirror of the cavity), we can achieve e ~1 and 
a maximum fiber coupling efficiency of over 35 %. A 
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FIG. 3. Probability of collecting photons with a fiber cavity 
into the mode of a single mode fiber as a function of cavity 
length (l). As the cavity length approaches the RoC of the 
spherical mirror, the beam waist at the fiber tip shrinks. Tf 
= 9,500 ppm and T^ + Lf + Le= 1500 ppm are assumed. The 
radius rion is smallest for a cavity length of 4.9995 mm while 
best mode overlap between cavity and fiber is achieved for a 
length of 4.9999 mm. 

comparison with a cavity constructed using two concave 
mirrors is presented in Appendix 1X1 

In summary, we expect to reach (g,K,7)/27r sa 
(18, 26, 10) MHz using a flber tip cavity. With pulsed 
excitation, one can achieve more than 30 % coupling effi- 
ciency into the single mode fiber if the minimum available 
RoC is about 5 mm, and the average rate of single photon 
emission process is on the order of max(l/(7, 1/k). 

C. Excitation by stimulated Raman adiabatic 
passage 

In an atomic system where the initial state \i) and the 
final states |/i), I/2) are different as shown in Fig. [l](b) 
and can be connected by a resonantly enhanced Raman 
transition, the STIRAP process can be used to generate 
a single photon in the cavity mode [l^, [l3| • 

For this process, a laser beam is resonant with one 
branch of the Raman transition (|i) <^ |e)) with Rabi fre- 
quency rip, and the cavity vacuum stimulates the emis- 
sion of the photon on the other branch (|/) |e)) with 
vacuum Rabi frequency 2g. As fip is slowly increased 
from zero, the ion adiabatically follows the dark state of 
the Hamiltonian starting from \i) to |/). Because the 
excited state |e) is not populated, spontaneous emission 
does not occur. The photon collection probability Pcavity 
is 2C/(2C+ 1). If used together with the fiber cavity ge- 
ometry discussed above, Pcavity ~ 0.56 and Pf ibcr ~ 0.40 
can be reached. 

If one is interested in an ideal single photon source, 
STIRAP is a better choice because of its larger photon 
collection efficiency per trial. Furthermore, in contrast to 
the fast excitation, the pulse shape of the emitted photon 
can be controlled by ^lp{t). For the protocol to gener- 
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ate entangled ion pair based on entanglement swapping 
between a photon pair, successful generation rate of the 
photon pair is the relevant figure of merit, and is pro- 
portional to Pfijj^j. ■ Rrep, where Rrep is repetition rate of 
excitation process. Therefore in addition to the availabil- 
ity of the necessary atomic states, the optimal excitation 
method depends on the details of the experimental setup, 
such as the time required to prepare the initial state as 
well as cavity QED parameters. 

To estimate the time required to complete a single STT 
RAP photon generation, we need to add the time neces- 
sary to prepare the initial state \i) as well as the time 
for the adiabatic evolution and/or emission. Adiabatic- 
ity requires that flp ^ g^, and if i7p increases linearly 
over time T, T ^ 1/g. In the strong coupling regime 
(g ^ K, 7), slow ramping of flp guarantees that the quan- 
tum state of the system will be adiabatically transferred 
from \i,n = 0) to \ f,n = 1) without spontaneous emis- 
sion where n is the number of photons in the cavity mode. 
The single photon then leaves the cavity with a rate of 
K. The Rabi rate ilp should remain at its maximum 
value until the single photon leaves the cavity. Therefore 
with strong coupling, the total photon generation time 
includes 0{1/k) as well as time necessary for adiabatic 
process (3> In the weak coupling regime, the cavity 

decay process will occur much faster than the adiabatic 
transfer, and therefore the minimum time is mainly lim- 
ited by T 1/g. Because Yb ion-cavity system lies 
in weak coupling regime, STIRAP process will require 
much more time to extract a photon compared to fast 
excitation case. In the rest of this paper, we will mainly 
consider fast excitation to compare the cavity result with 
high NA optics in the following section. 



III. COLLECTION OF PHOTONS WITH HIGH 
NUMERICAL APERTURE OPTICS 



As an alternative method to collect light from a single 
trapped ion, we consider the use of high NA collection 
optics whose geometric aberrations are well-corrected to 
achieve diffraction limited coupling performance into a 
single mode fiber. Reflective optics is frequently used as 
an alternative to refractive optics to realize a large NA 
[30l - l36j . A parabolic mirror can perfectly collimate light 
emitted from a point source, and has been studied ex- 
tensively [3^-l36^. In this section we will investigate the 
properties of light reflected by a spherical mirror and nu- 
merically calculate the collection probability into a sin- 
gle mode fiber. Although the analysis is presented for 
refiective optics, the general conclusions apply to other 
approaches if similar NA and aberration control can be 
realized. 




FIG. 4. Definitions of coordinate systems and symbols with 
respect to spherical mirror, (a) Ion is located at the origin 
of coordinate system, and polar (cylindrical) coordinates are 
used to describe the optical field before (after) refiection. 
is the polar angle with respect to S-axis, and is azimuthal 
angle from i-axis. (b) 6max is the maximum half opening 
angle corresponding to the mirror aperture, /i is a distance to 
the analysis plane from the center of sphere. Surface ion trap 
can be fabricated on the substrate area surronding mirror. 



A. Dipole emission and Gaussian mode 

Distinguishing light from tt and a transitions is essen- 
tial for achieving high fidelities in remote entanglement 
protocols. Therefore, the modes of photons from tt and 
a transitions have to be orthogonal after collection. For 
large numerical aperture, this is only possible if we choose 
B parallel to the optical axis z. The optical fields of pho- 
tons emitted by the ion via the three possible transitions 
are given by [35| 

£^A;=i,Am=o = \/-^ sin 6*6', 

r V OTT 

^A/=i,Am=±i = '-^e^'\^-{±cose6 + i^). (6) 
r V lOTT 

Once it is reflected by a curved surface, spherical waves 
from an ion become paraxial waves, and the reflected 
optical field will be described in cylindrical coordinates 
shown in Figure |4] (a) . If the curved surface is parabolic 
and the trapped ion is located at the focus of paraboloid, 
all the reflected rays propagate parallel to the z-axis and 
the traveling wavefronts are flat. 

For a spherical mirror, most of the reflected rays prop- 
agate with radially varying phase delay and non-zero 
angle with respect to the z-axis due to spherical aber- 
ration, even if the ion is trapped at the nominal focus 
(midpoint between the center of sphere and the vertex 
of the spherical surface). In this case, rays cross each 
other and the optical field distribution in p-coordinate 
changes as it propagates along the z-direction [Figure |4] 
(b)]. The probability of collecting a photon into a single 
mode fiber can be estimated by numerically calculating 
the mode matching factor [Eq. ([S])] between the reflected 
optical fleld distribution (Ri^m) incident on the flber and 
the Gaussian flber mode {G fiber) by setting E = Ri.m- 
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FIG. 5. Maximum collection probability of photons emitted 
by a CT-transition (AZ = l,Am = ±1) and a 7r-transition 
(Ai = l,Am = 0). For the yr-transition calculation, we as- 
sume that phase e"^ is added to the collimated field. We also 
plotted the amount of light collimated by parabolic mirror 
before being collected by single mode fiber, po is the radius 
of mirror aperture and / is the focal length of the parabolic 
mirror. Inset: Spherical probability distribution of photons 
emitted by different transitions with respect to the quanti- 
zation axis. Distance from the origin is proportional to the 
probability of photons propagating in the corresponding an- 
gle. 



The collection probabilities of light from a single ion 
onto a large area detector (referred to as "collimated 
case" ) and into a single mode fiber ( "fiber coupled case" ) 
are shown in Figure [5] as a function of the mirror aper- 




FIG. 6. Matching collimated ray distribution to Gaussian 
beam with finite wavefront radius of curvature. For clarity, 
figures are not drawn to scale. 



Luo et al. ,35j analyzed the case of a parabolic mir- 
ror and calculated the overlap of the reflected mode 
with a Gaussian mode. They showed that the collec- 
tion probability of a photon emitted by the 7r-transition 
{Al — l,Am — 0) vanishes due to the rotational cylin- 
drical symmetry of polarization, while the collection effi- 
ciency for (T-transition {Al — 1, Am — ±1) approaches 50 
% as the mirror aperture (po) becomes much larger than 
the focal length (/) of the parabolic mirror. Provided the 
axis of single mode fiber is collinear with the quantiza- 
tion axis, a single mode fiber can be used as transition- 
selective filter which maps photons generated by and 
cr~-transitions into two orthogonal circular polarizations 
{{x =F ) respectively, and rejects photons gener- 

ated by the tt transition. For the frequency qubit shown 
in Fig. [l](d), it would be interesting to collect photons 
from TT-transitions and reject a transition photons [l^ . 
We can implement this orthogonal filter by adding an ad- 
ditional (/)-dependent phase (e*"^) to the optical field after 
it is reflected by the mirror [Eq. (33) and (35) of Ref. 
[35}]. A phase that depends on the azimuthal angle can 
be added using computer generated holograms [s^, HI] , 
diffractive optics 39, 40] phase plates |41„[42], or spa- 
tial light modulators j^J, Sj] . Therefore we assume that 
whenever we intend to collect light from 7r-transitions, 
an additional phase factor e*"^ is added to the (near-) 
collimated field. 



ture radius (po) for both a- and 7r-transitions. When 
Pq is smaller than /, the fiber coupling efficiency of 
(T-transition photons is much larger than that of tt- 
transition photons due to the dipole emission pattern 
[inset in Figure [5]. When po is much larger than /, spa- 
tial overlap with fiber mode becomes less sensitive to the 
details of the emission pattern but the spatially varying 
polarization distribution limits the overall coupling effi- 
ciencies of both transitions near 50 %. 

The integration of an ion trap with large aperture 
parabolic mirrors has been considered [sol] , In this work, 
we will concentrate our analysis on the collection proba- 
bility one can achieve with a spherical mirror, which can 
be integrated with a surface trap, as shown schematically 
in Figure|3](b) [3l|. As spherical mirrors feature substan- 
tial aberrations at large NA, this analysis also provides 
strategies for effectively managing aberrations present in 
most practical collection optics. When the ion is trapped 
at the nominal focus of the spherical mirror, the rays col- 
limated by a spherical mirror are not perfectly parallel 
to the optical axis but converge initially and eventually 
diverge after crossing the optical axis (Figure H]). The 
resulting far-field wavefronts are close to spherical with 
residual optical path difference (OPD) as shown in Fig. 
ini To estimate the probability of coupling this mode to 
an optical fiber, we calculate the mode overlap to a Gaus- 
sian mode with spherical wavefronts: 



= — J-cxp(iV')exp ( - 

K7 V TT \ W''(Z) 



exp 



\2R{z) 



(ax + Py), 



(7) 



where R{z) is the radius of curvarture (RoC) of the wave- 
front, is the beam radius at the measurement plane, 
and "0 = arctan(— z/2z7?) is Gouy phase. The residual 
OPD between the wavefront of collimated light field and 
the spherical wavefronts of the Gaussian mode limits the 
maximum coupling efficiency one can achieve with this 
method. The coupling efficiency can be recovered to val- 
ues similar to the ideal parabolic case if the residual OPD 
is reduced to below the diffraction limit. 



B. Design examples 

The residual OPD (expressed in units of wavelength A 
= 369.5 nm) of the rays refiecting off a micromirror with 
RoC R = 160 /Ltm, measured at a location /i = 50 mm 
compared to a best-matching Gaussian mode is plotted 
in Figure [T] The ion is assumed to be located at the 
focal point of the spherical mirror. The two horizontal 
gray dashed lines show the limits of Rayleigh criterion 
representing the diffraction limit (OPD < ±A/4), below 
which residual OPD does not lead to reduced coupling 
efficiency The black dashed line shows the resid- 

ual OPD for a spherical micromirror with opening angle 
{Omax) of 32°, where 10.5 % of tr-transition light falls 
within the aperture and is collimated by the mirror. In 
this case, the OPD remains within the Rayleigh crite- 
rion, and a numerical calculation estimates that 6.2 % 
of photons from cr-transitions can be matched to a single 
mode fiber. In contrast, the light field generated by the 
TT-transition has little overlap with a Gaussian mode for 
Omax — 32° due to the emission pattern. 

When Omax is extended to 48° to increase the amount 
of collimated light from 7r-transition to 7.3 %, the maxi- 
mum OPD becomes larger than A due to spherical aberra- 
tions (dotted line in Figure!?]) and the maximum coupling 
efficiency is still limited to 0.4 %. The residual OPD must 
be compensated in order to recover the fiber coupling ef- 
ficiency. If the distribution of the OPD is known, one can 
utilize a phase plate [4l|, HI] or a spatial light modulator 
(43I li^ l to introduce compensating OPD to cancel out 
the residual OPD. Inset of Figure [7] shows an example 
of a phase plate which can be fabricated on a 500 /im- 
thick fused silica substrate to compensate the residual 
OPD for the 7r-transition case. Using this phase correc- 
tor plate, the coupling efficiency into the fiber mode can 
be enhanced to 3.0 %. We used /i = 50 mm as a typical 
location to place phase correcting optics, but the general 
conclusion is independent of where the mode matching is 
calculated. We also explored optimization of ion height 
along the optical axis, and concluded that no substantial 
improvement can be achieved. 

Figure |S] shows the probability of coupling emitted pho- 
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FIG. 7. Example of residual OPD expressed in units of wave- 
length (A = 369.5 nm). The RoC of matching wavefront is 
chosen to maximize the mode overlap rather than minimizing 
the maximum OPD. According to Rayleigh wavefront crite- 
ria [4^, maximum OPD should be less than A/4 to avoid 
the effect of non-zero OPD and this range is marked by two 
dashed lines. The dark dashed hne and dotted line shows the 
residual OPD (measured at a location /i = 50 mm) for the 
rays reflected by a spherical mirror {R — 160 fim) with Omax 
— 32° and Omax ~ 48°, respectively. In the inset, thickness 
variation of phase plate is plotted, designed to compensate 
OPD shown by the dotted line. 500pim-thick fused silica sub- 
strate is assumed in this design. With this phase correcting 
plate placed at ft = 50 mm, residual OPD at ft = 50 mm 
is dramatically suppressed. Solid line is OPD of phase cor- 
rected wavefront measured at ft = 100 mm, which shows that 
OPD barely changes for additional propagation once OPD is 
compensated. 



tons into a single mode fiber using a spherical mirror and 
an optimized phase corrector plate for each value of 0max ■ 
With optimum phase correction, the collection probabil- 
ity as a function of the opening angle O^ax of a spherical 
mirror grows similar to that of a parabolic mirror until 
Omax reaches ^ 50° where matching of slightly diverg- 
ing rays to a Gaussian mode starts to fail. Beyond this 
limit, we need to correct the angle of collimated rays in 
addition to the phase which is beyond the scope of this 
paper. Likewise, 7r-transition coupling efficiency could 
be improved to over 3.0 % from 0.4 %, approaching the 
parabolic mirror limit if all aberrations are properly com- 
pensated. 

Use of a microfabricated mirror integrated under a 
surface trap has been proposed for efficient collection of 
emitted photons for qubit state detection ^ . In that ap- 
proach, however, the aberrations of the microfabricated 
mirror cannot be controlled well enough to produce a fiat 
wavefront. One important advantage of using a microfab- 
ricated mirror compared to a conventional macroscopic 
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FIG. 8. Fiber coupling efficiency of a photon emitted by a- 
and TT-transition using a spherical mirror with a phase correc- 
tor plate when R — 160/^m. We also plotted fiber coupling 
efficiency by parabolic mirror for comparison. This plot shows 
that most of aberration of spherical mirror can be corrected 
by phase plate when 6max < 50°. 



Mirror RoC (/im) 


1600 


160 


16 


Transition 


a 


TV 


a 


TT 


0" 


TV 


Collimated light by mirror (%) 
Fiber coupling probability (%) 


21.2 
<0.1 


7.3 
<0.1 


21.2 
6.2 


7.3 

0.4 


21.2 
15.8 


7.3 
2.8 



TABLE I. Estimation of fiber coupling efficiency of a and 
TV transition light by ideal spherical mirrors of different di- 
mensions. 6max is set to 48°. In our simulation, we added 
an ideal hyperbolic lens in the beam path, which can trans- 
form the matched Gaussian mode into UV fiber mode {wo ^ 
1.54 fim) and performed full ray tracing from the ion position 
to the fiber plane to evaluate the mode overlap integral. 



mirror [33| is that the geometric aberrations scale as the 
dimensions of the optics are reduced [i^. The residual 
OPD is reduced while the wavelength (A) remains the 
same eventually achieving diffraction limit. Table |T] sum- 
marizes photon coupling efficiency using a single spher- 
ical mirror with the same large NA {9max ~ 48°) with- 
out any additional aberration correction. Last column 
for RoC = 16 fj,m shows that there will be a significant 
improvement if the dimensions of the mirror and the po- 
sition of the ion are shrunk by a factor of 10. For such 
micro optics, the detailed shape of the reflector has little 
impact on the collection efficiency. On the other hand, 
if the dimensions of mirror are increased by a factor of 
10, the collection probability will vanish unless there is 
proper phase compensating optics (first column) . There- 
fore the microfabricated mirror with RoC ^ 100 fj,m can 
increase the fiber collection efficiency drastically. How- 
ever due to technical limitations of the ion height, it is 
impractical to use mirrors with RoC much smaller than 
~ 100 fim at present day. 



C. Alternative approaches 

The analysis in this section shows that efficient cou- 
pling of light into a single mode fiber can be achieved 
provided one can correct the aberration of a large NA col- 
lection optics to conform to the wavefront of a Gaussian 
beam, to below the diffraction limit. Practical means to 
realize such high NA optics with well-corrected aberra- 
tion include (1) design of multi-element refractive lenses, 
(2) fabrication of paraboloidal mirrors (using e.g. dia- 
mond turning technology) , (3) spherical mirror with cor- 
rector plates to compensate for the residual aberration, 
and (4) Fresnel refiectors or diffractive optical compo- 
nents. 

Fresnel lenses have been investigated as a way to collect 
light with large NA without spherical aberration [13, 13 
because it adds the necessary phase shift to a continu- 
ous wavefront. Using the similar design principle, we can 
consider integrating Fresnel refiector with an ion trap in 
place of the micro mirror shown in Figured (b). One 
critical consideration for Fresnel diffractive optical ele- 
ments is the diffraction efficiency that depends on the 
grating periods [47] . The Fresnel refiector deserves more 
careful analysis for use in aberratioin-free fiber coupling 
applications. 



IV. APPLICATION TO REMOTE 
ENTANGLEMENT GENERATION 

Based on our estimation of photon coupling efficiency 
into a single mode fiber, we can calculate the expected en- 
tangled ion pair generation rate considering practical ex- 
perimental constraints. Table |TT] compares the expected 
entanglement generation rate of various approaches with 
existing experimental results. Matsukevich et al. (isj 
generated entangled ion pairs by using polarization qubit 
[Figure [T] (b)]. Because state preparation step requires 
only optical pumping into ^S'i/2|0,0) without any qubit 
rotation [l^, the same procedure can be repeated at 520 
kHz repetition rate, and the entangled pair generation 
rate is analyzed in the first column of Table [Hi When 
the remote entanglement is established by coincident de- 
tection of a photon pair, the ionic qubit is encoded in the 
Zeeman substates ^51/211,^1) which are susceptible to 
magnetic field fluctuation and therefore these two states 
should be transferred to atomic clock states ^51/2! ± 1, 0) 
by microwave transitions or Raman transitions. The pho- 
tonic polarization qubit is also susceptible to decoherence 
from birefringence in the vacuum window, downstream 
optics and single mode fibers. Therefore this method 
requires careful control of potential decoherence mech- 
anisms to achieve high fidelity. Olmschenk et al. [l3| 
used the first-order magnetic field-insensitive hyperfine 
clock state qubit accompanied by tt transition as shown 
in Figure [T] (d), and collected emitted photons from a di- 
rection orthogonal to the quantization axis to maximize 
collection of photons generated by tt transition [inset of 
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Figure [5] . To generate ion-photon entangled state, the 
ion should be prepared in a superposition state before ex- 
citation, which requires rotation of clock states using mi- 
crowave transition and limits the overall repetition rate 
to 75 kHz. The analysis of this experiment is listed in the 
second column of Table HIl The limitation of this collec- 
tion geometry is the increased overlap between different 
polarizations from a and tt transitions, which reduces fi- 
delity as NA is increased [l4| . 

From the third to the fifth column in Table |TT] shows 
the expected entangled ion pair generation rates when 
the numerical aperture is increased to 0.6 using either 
refiective optics or refractive optics discussed in Section 
mil We considered several practical improvements that 
can be made over the existing experiments. Although 
previous experiments [T2| - [l^ detected only the singlet 
state out of four Bell states, one of triplet states can also 
be distinguished from the other two. Therefore we as- 
sume that probability of Bell state measurement is 1/2 
rather than 1/4. In previous experiments, various losses 
in photon collection efficiency caused by non-ideal opti- 
cal components further reduced the overall efficiency by 
a factor of ^0.25. We estimate that we can use better op- 
tical components to increase this efficiency to more than 
60 % based on realistic analysis. Third column analyzes 
experimental setup similar to the first column [l3| when 
the numerical aperture is increased to 0.6. This setup 
can achieve the highest expected generation rate among 
the three cases, but requires decoherence control simi- 
lar to the first column setup. Forth column lists values 
we can achieve using similar setup to the second column 
[3| , but we assume that optic axis of collection optics is 
parallel to the quantization axis compared to orthogonal 
setup used in the second column. As was discussed in 
section IIIIl if a phase plate is used in the parallel con- 
figuration, only the tt transition can be coupled into a 
single mode fiber and the fidelity will remain high even 
with very large numerical aperture. However, NA of 0.6 
is not large enough to show significant improvement due 
to non-ideal emission pattern of 7r-transition in this col- 
lection direction. Fifth column estimates the generation 
rate which can be achieved with a scheme in Figure [T](c) 
[l^ . This approach has advantages of high fidelity due 
to frequency qubit and a transition which allows high 
collection efficiency. However, to achieve near unity ex- 
citation probability, the initial state should be prepared 
to either a superposition of the clock states or one of 
^'5'i/2|lj which limits the repetition rate around 75 
kHz. 

The last column shows the entanglement generation 
rate which can be achieved by an optical cavity discussed 
in Section ini In schemes using optical cavities, the FSR 
of cavity should be matched to cither 12.6 GHz or 14.7 
GHz to collect photons encoding frequency qubits. These 
constraints lead to minimum cavity length of ^10 mm, 
precision control of RoC of the spherical mirror down 
to /xm level, and an upper limit on the cavity finesse of 
~ 700 to cover 19.6 MHz hncwidth of the '^Pxji state. 



Therefore frequency qubit is not compatible with optical 
cavity schemes. However, we can achieve highest entan- 
glement generation rate compared to other approaches 
using polarization qubit or time-bin qubit. In this case, 
the remaining challenge will be the minimization of de- 
coherence similar to all the other schemes involving the 
polarization qubit. 

Based on the above analysis, collection of photons gen- 
erated by <7-transition along quantization axis is the op- 
timal choice for both the large NA system and the cavity 
system. Cavity system prefers polarization qubit as the 
photonic qubit of choice, but it can achieve the highest 
entanglement generation rate. For large NA system, fre- 
quency qubit has lower generation rate than polarization 
qubit mainly due to slow initial state preparation. In 
principle, this limitation can be overcome with fast mi- 
crowave transition or Raman transition. A custom lens 
system designed for optimal chamber configuration can 
achieve a NA of more than 0.6 with diffraction-limited 
performance. Even though the maximum collection prob- 
ability with free space optics is smaller, free space optics 
has the advantage that it can be realized with off-the- 
shelf technology today. Therefore, a-transition-based fre- 
quency qubit system with large numerical aperture lens 
can be a good compromise between high fidelity and high 
throughput. 



V. CONCLUSIONS 

In this paper, we discussed various ways to enhance 
collection of photons emitted by a single atomic ion into 
a single-mode optical fiber. Based on this analysis, we 
conclude that the generation rate of entangled ion pairs 
through entanglement swapping of photons can be dra- 
matically enhanced beyond the current values using ei- 
ther a high NA collection optics or a carefully designed 
optical cavity. Such enhancement is expected to enable 
the realization of scalable quantum multi-processor and 
quantum repeaters based on nodes of trapped ions con- 
nected via photonic network. 

This research was funded by Office of the Director of 
National Intelligence and Intelligence Advanced Research 
Projects Activity through Army Research Ofl[ice. 

Appendix A: Cavity with two concave mirrors 

In this Appendix, we estimate the maximum collec- 
tion efficiency one can achieve with a cavity made of 
two curved mirrors with an ion trapped in the center. 
Compared to the fiber cavity, ion can be trapped at 
the beam waist of the cavity mode, and we can esti- 
mate the lower bound and the upper bound of Pcavity 
In the near-concentric limit, we estimate the optimal 
Wo to be < uiQ < \/ A/tt-^/RoC • A/2, which corre- 
sponds to 1 > 1/(1 -I- {wo/rc)"^) > 0.87. Combined with 
1/(1 -|-2RoC/Zc) w 0.53 which is possible with symmetric 
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Numerical aperture 


0.23 


0.6 


cavity 


Optic a^xis vs. B-field 


II 


_L 


II 


II 


II 


II 


Type of atomic transition 


a 


TT 


o 


TT 


a 


a 


Photonic qubit 


Pol. 


Preq. 


Pol. 


Preq. 


Preq. 


Pol. 


Detector efficiency 


0.2 


0.2 


0.3 


0.3 


0.3 


0.3 


Decay via right transition 


2/3 


1/3 


2/3 


1/3 


2/3 


n.a. 


Collected fraction 


0.0198 


0.0198 


0.136 


0.028 


0.136 


0.337 


Mode overlap 


0.82 


0.82 


0.85 


0.32 


0.85 


0.95 


Misc. efficiency 


0.26 


0.25 


0.63 


0.62 


0.63 


0.63 


Single Photon collection efficiency (%) 


0.0554 


0.0269 


1.46 


0.0547 


1.46 


6.04 


Bell state identification 


0.25 


0.25 


0.5 


0.5 


0.5 


0.5 


Coincidence efficiency 


7.7 X 10"* 


1.8 X 10"* 


1.1 X 10"" 


1.5 X lO"'^ 


1.1 X 10"* 


1.8 X 10"^ 


Repetition rate (kHz) 


520 


75 


500 


75 


75 


500 


Entanglement rate (Hz) 


0.04 


0.0014 


53 


0.011 


8 


913 



TABLE H. Comparison of entanglement generation rate based on different setup. Theoretically maximum two types of Bell 
states encoded in a photon pair can be unambiguously identified from the others, and we include this identification efficiency 
on the 5th row. When atom decays to the ground state, some of the decay channels are undesircd and we included the fraction 
of desired decay channel in 7th row. Collected fraction corresponds to the fraction of collimatcd light when the photous are 
collected with reflective or refractive optics and r* • PcavUy for the case of the cavity system. Experimental condition of each 
column is described in the main text. 



geometry, we can estimate that 0.53 > PcavUy > 0.46 ity. Mode matching has to be performed to couple these 
which is comparable to Pcavity ~ 0.41 of the fiber cav- photons into a single mode fiber, and it is not easy to 

integrate this geometry with surface traps. 
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